Reconstructing the density operator via simple projectors 
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We describe the representation of arbitrary density operators in terms of expectation values of simple projection operators. 
Two representations are presented which yield non-recursive schemes for experimentally determining the density operator of 
any quantum system. We suggest a possible experimental implementation in quantum optics. 
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I. INTRODUCTION 

In the realm of quantum theory a state of a physical 
system is most generally expressed by its density operator 
Q. Knowledge of this operator gives complete informa- 
tion of the quantum state. Schemes have been proposed 
in a number of fields in quantum physics to determine g 
from measurements either explicitely 5 , |^ | or in- 

directly via quasiprobability distributions ME,^^ for 
mixed states and also for pure states only |1^ 1^ . 

In this Letter we describe a general method of repre- 
senting any density operator g in terms of expectation 
values of simple projection operators. Since the expecta- 
tion values of projectors can, in principle, be determined 
experimentally this approach leads to schemes for exper- 
imentally determining the density operator. 

Our approach differs from previously proposed schemes 
in quantum optics for determining the density operator 
in its use of simple projectors which project onto a sin- 
gle or a linear superposition of two basis states [§,1). 
We place an emphasis on a 'minimalistic' representation 
which comprises the least number of projection opera- 
tors and thus leads to the most efficient scheme, it is a 
generalization of the previous considerations in |^,^. 

The plan of the paper is as follows. First we intro- 
duce the general idea of our approach, then, in section 
II, we cast it into two specific representations and de- 
scribe their relative virtues. In section III we describe a 
quantum optical implementation and we end with a dis- 
cussion in section IV. 

Let us assume that the Hilbert space representing the 
states of the physical system is of countable dimension N 
and let |m) for m = 1, A'^ be any conveniently chosen 
orthonormal basis of the space. In cases where the space 
is infinite in dimension, all expressions containing N here 
and in the following are infinite also. Our primary aim is 
to represent the iV^ — 1 independent density matrix ele- 
ments gnm = {n\g\m) in terms of the expectation values 



of simple projection operators. Clearly the matrix ele- 
ments cannot be expressed solely in terms of the TV — 1 
independent expectation values gmm = {\fn){'m\) of the 
set of the N base state projectors \m){m\ because the vi- 
tal phase information of the coherences, i.e. the complex 
nature of the off-diagonal elements gnm for n ^ m, can 
not be derived from the diagonal elements alone. 

The simplest possible generalization of the base state 
projectors is given by the set of projection operators 
which project onto linear combinations of two base states, 
e.g. onto ci\n) -|-C2|to). The expectation value of such 
projectors represent the phase information of the coher- 
ences in its most elementary form. We show that one can 
express g in terms of expectation values of such projectors 
and how to implement it quantum optically. 



II. TWO REPRESENTATIONS 

For simplicity let us consider the {\n), |m)}-subspace 
which is spanned by any two basis vectors \n) and |m) 
for n ^ m and define the state 



n)+a\m)), (1) 
is a normalization constant 



where Na = l/y/TT\a, 
and a = |a|e'" is a nonzero coefficient. A corresponding 
nomenclature is used for a second, different state of the 
same subspace \h) = Ab(|n) -|-6|m)), where b = \h\e^^ ^ a. 
We defer making any further restriction on the values of 
a and 6, to guarantee independence of the expectation 
values of the corresponding projectors 



A = \a){a\ , 
B^\b){b\ , 



(2) 



until later. 

Let us assume the measurements yielding the expec- 
tation values of the projectors \m){m\, A and B 
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have been performed ||T^ . The first two expectation val- 
ues are simply the diagonal elements p„„, gmm- We can 
combine these expectation values conveniently as 



M|„) = Ti{qA} - iVf (f?„„ + \a\^g,nm) 



(3) 



where Tr is the trace operation and -/^^|a) stands for the 
result associated with a measurement of the projector 
A. A corresponding expression is obtained for the result 
M|5) associated with the projector B . Let us write Qnm 
in terms of its real and imaginary parts Qnm ^ R + iJ 
and let us define 



ia) 



2\a\m 



— R cos a — J sin a 



and mih) = ^j^j-^ = i?cos/3 — Jsin/3 
Solving these equations for R and J yields 



(4) 



1 / sin j3 — sin a 
sin(^-a) \cos/3 -cos a 



'\a) 



(5) 



Clearly this requires [3 — a ^ kn, where k is any integer. 
This gives the only restriction on the values of a and b 
aside from the trivial requirement that a ^ ^ b. Ap- 
plying the outlined procedure to the {\n), |m)}-subspaces 
for 1 < n < m < N allows us to represent g in terms of 
expectation values of N'^ ~ 1 different projectors, due to 
the condition Tr g = 1. Note that this scheme is intrin- 
sically non-recursive. 

We call this the 'minimal' representation as it requires 
this least possible number of projection operators to rep- 
resent a general density operator and also because it puts 
almost no restrictions on the states forming the projec- 
tors, namely on the coefficients a and b of Eq. 

Though mathematically satisfactory the minimal rep- 
resentation would be sensitive to experimental errors in 
a physical implementation. This sensitivity however is 
minimized using sensitivity optimized states, i.e. choos- 
ing \a\ = |6| = 1 and b = ± i a This sensitivity can 
be further reduced employing three or more (redundant) 
states. Let us for example look at the case of one more 
projector state |c) = Nc{\n) + c|m)) where c = |c|e'^ in 
each {n, m}-subspace. We find that 



'\c) 



m|a) sin(/3 - 7) - m|6) sin(Q; - 7) 
sin(/3 — a) 



(6) 



where m^^.^ is given by Eqs. (Q,^ with \b) replaced with 
|c). Provided the differences a — 7, /3 — 7 and a — (3 be- 
tween the phase angles of the states \a), \b) and |c) are 
not multiples of tt the overparameterization introduced 



by the extra state can be used to reduce the effect of ex- 
perimental errors. For example, one could estimate true 
values of m\a) , ™|&) and m|c) as the point (x, y, z) on the 
surface z{x, y) = [x sin(/3 — 7) — y sin(a — 7)]/ sin(/3 — a) 
which is closest to the point {x,y,z) where x, y and 
z are the experimentally measured values of mi^^^, m^i,'^ 
and m|c) . 

One may still go one step further and consider the 
particular quadruplet of states 

|ar>^^(l^)±H), 



\br)^^{\n)±i\m)), 



(7) 



for n,m = l,2,...,N. We mention in passing that 
all such states are normalized except for n = m for 
which la!;:") = V2\n) and |a™),|6'i") = 0. The set 
{la^-™), : m,n — 1,...,N} is an overcomplete ba- 

sis of the Hilbert space. Let the projection operators |T^ ] 
which project onto these states be A™ = |a5™)(a5™|, 
BJ™=|65^™)(6!^"'|, defined in analogy to Eq. (|). The 
expectation values of the 2N'^ — N different projectors 
i P?! for n < m suffice to represent an arbitrary matrix 
element of g as 



KB', 



B!!™)]}, (8) 



a form that has already been derived in f^^- Now the 
projectors can be combined to form operators _R"™, J"™ 
defined as 

jnm ^ (^„™ _ B'l"')/V2 = i{\n){m\ - \m)(n\)/V2 , (9) 
fulfilling the orthogonality relations 

Tr{i?"™7??'n = {Sn,pSm,g + Snjm,p) , 

and Tr{i?"™ = , (10) 

for n,m,p,q = 1,...,N, where S„^rn is the Kronecker 
delta. The set {i?™", J""" : n < m} constitutes a 
complete basis set of operators. This operator basis 
gives an unique expansion of any operator Q as 



N m-1 



q = v2{y.y: 



m— 2 n—1 



Jnm 'J 



(11) 



with r^m = Tr{QR"'''}/V2 = {Qmn + Qnm)/'2 and 

Jnm = Tr{QJ™"}/V2 = [Qmn " Qnm)i/2. If Q is a 
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hermitian operator r„m and jnm are the real and imagi- 
nary parts of the matrix elements Qnm = {nlQlm). 

Fano introduced the idea of expanding the density ma- 
trix in terms of an orthogonal operator basis Q, hence 
we call this an 'operator basis' representation. We intro- 
duced this representation for its mathematical properties 
rather than its physical contents. Let us note that the 
sensitivity optimized states mentioned before Eq.(||) can 
analogously be cast into this kind of orthogonal operator 
basis, in this sense the operator basis representation is 
contained in the minimal one. 



III. QUANTUM OPTICAL REALIZATION 

Next we describe a possible experimental scheme for 
the reconstruction of a density operator describing the 
state of a single optical field mode ||l8| . It is a straightfor- 
ward matter to generalize this to several optical modes. 
We use the Fock state basis in which the numbers of pho- 
tons in the mode under consideration label the states 
{|m) : m = 0,1,2,...}. Our task is to show that the 
expectation values of the corresponding projection op- 
erators AjJ^, A, etc. can be obtained experimentally. 
We note from the outset that the experimentally difhcult 
part of the scheme at present is the preparation of coher- 
ent superpositions of two Fock states. However, in the 
light of recent theoretical ||l^,^ and experimental results 
pl[ , it is clear that the problem of the preparation of the 
probe field can and will be solved. 

Thus, since this is not a fundamental difficulty we as- 
sume in the following that such superposition states are 
available. 

The expectation value of the projection operators in 
the representations can be determined using the experi- 
mental setup depicted in Fig. 1 as follows. 



If you are interested in the figure, please contact me via 
email or fax and send your fax number, then I'll send the 
picture by fax. 
Ole Steuernagel 

email : ole@photon.fta-berlin.de 
fax: ++49 (germany) -30- 6392 3990 



FIG. 1. The setup of our proposed quantum optical 
scheme. Light from a common field source is fed into a device 
generating the probe field \a){a\ and a device that generates 
the signal field q. The probe and signal fields, which are la- 
belled 1 and 2, respectively, are then entangled at the last 
beam splitter and analyzed by the photodetectors / and //. 
The use of a common source ensures that the probe and signal 
fields oscillate at the same frequency. 



A probe field is prepared in a particular state and 
fed into port 1 of the beam splitter, the signal field pre- 
pared in the (unknown) state g is fed into port 2. The 
joint photon number probability distribution of the out- 
put ports of the beam splitter is obtained from the pho- 
toelectron statistics produced in the photodetectors / 
and // for many repetitions of the experiment, let us 
note that multiphoton coincidence counts together with 
quantum efficiencies above 70% have been demonstrated 
experimentally |^,^ . If one chooses a method that de- 
tects single photons with more than 50% quantum effi- 
ciency the photon number probability distribution can 
be recovered from the measurements using the inverse 
Bernoulli transformation discussed by Lee |24[| . 

Furthermore a new method developed by Munroe et 
al. |2^ allows to measure the photon-number statis- 
tics from the phase-averaged quadrature-field distribu- 
tion with single photon and ultrahigh time resolution of 
the order of 300 fs. Employing the corresponding recon- 
struction schemes |^ this method yields almost perfect 
photon number statistics. 

Hence we may restrict our considerations to the 'true' 
joint photon probability distribution {p, q) for p and 
q photons measured by (ideal) photodetectors / and //, 
respectively, which is given by 



p+q p+q 



m 



n'—Q m'—Q 



xAp{n',p + q — n')A*{m',p + q — m') . 



(12) 



Here Ap(h', fj.) represents the probability amplitude of 
finding mode / in the Fock state \p)i if modes 1 and 
2 are in the product Fock state \v)i\fi)2 and is given by 



,i<,3p(p-i/) 



A>,^) = (-l)^/fcV^e-(--) e>^ 

X E E(-l)' il) (fj ^r^^+^-'p'^-^+i 4+,p , (13) 
fe=0 (=0 \ / \ / 

where r, p are the transmittance and reflectance and Lpr, 
Lpp are the corresponding phase factors generated by the 
beam splitter as defined by Campos et al. in Ref. [ p7| . 
Inserting for the special probe field states |a"™) with 
n > m, see E q. (p] ), and relabeling p+g = N + n = M + m 
changes Eq. (|12|) into 



P|,„,„)(p,iV+n-p) = C(a'|e|a'), 



(14) 



where 



|a')- C-i/2 [A;{N,n)Na\N)+A;{N,n)Naa*\M)] 
and C= \Ap{N,n)Na\^ + \Ap{N,n)Naa\^ . 

(15) 

Again we assume that the diagonal elements are known, 
for example by the comparatively simple measurement of 
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the photo count distribution of the field alone. The same 
is assumed to be true for a, which is known from the state 
preparation process, we can thus, equivalently to Eq. (||), 
use P|Qnm^ (p, N + n — p) to determine a quantity 

M|,„„.)(iV,p) = 2Re{a qmn Ap{M,m)A;iN,Ti)} , (16) 

where 'Re' signifies the real part. Using a second linearly 
independent probe state jfe"™), by a procedure analo- 
gous to Eqns. and we obtain qmn- Thus we have 
translated the minimal representation into an experimen- 
tal scheme in quantum optics for determining the quan- 
tum state of light; the translation of the other represen- 
tations along similar lines is straightforward. 

It is interesting to note that the value of p in Eq. (|l^ ) 
can be chosen arbitrarily from the interval (0 < p < 
N + n). This gives A^-|-n-t-l different ways of determining 
the value of the quantity M\g^n.mj(N,p) in Eq. (|lq ). Also, 
since we require n — m — M — N in Eq. (|lj) the set 
of matrix elements g(k+N-M) k for A: = 0, 1, 2, ... can be 
determined from just two probability distributions P^a^m) 
and P|bnm^ for fixed values of n and m. And finally, since 
it is the difference n — m only that decides which set of 
matrix elements are determined this implementation is 
also redundant in the sense that the probe states |a**) 
with s ^ t + n — m are equivalent for t = 0, 1, 2, .... 

This scheme will give as many matrix elements of the 
density operator as desired and is limited only by exper- 
imental error and the ability to prepare the probe field 
in suitable two-Fock state superpositions. 

IV. DISCUSSION 

We examined the requirements for representing any 
density operator in terms of expectation values of simple 
projection operators. We gave two different represen- 
tations: the minimal representation which requires the 
least number of projectors and the operator basis repre- 
sentation which gives the expansion of any operator in 
terms of an operator basis. Our results are applicable 
to any physical system whose state space is of countable 
dimension TV which need not be finite. 

We showed how the expectation values could be de- 
termined experimentally for the case of a single mode of 
an optical field. An important point about our method is 
that it is not recursive, in contrast to some other methods 
for determining the density operator of the optical field 
P, p^p^ for which the calculation of all but a select few 
matrix elements involves the previously calculated values 
of other matrix elements and results in the accumulation 
of experimental errors. 
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